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Abstract 

We propose the time-dependent generalization of an 'ordinary' autonomous human 
biomechanics, in which total mechanical + biochemical energy is not conserved. We in- 
troduce a general framework for time-dependent biomechanics in terms of jet manifolds 
derived from the extended musculo-skeletal configuration manifold. The corresponding 
Riemannian geometrical evolution follows the Ricci flow diffusion. In particular, we 
show that the exponential-like decay of total biomechanical energy (due to exhaustion 
of biochemical resources) is closely related to the Ricci flowQ on the biomechanical con- 
figuration manifold. 

Keywords: Time-dependent biomechanics, extended configuration manifold, configu- 
ration bundle, jet manifolds, Ricci flow diffusion 

1 Introduction 

It is a well-known fact that most of Hamiltonian and Lagrangian dynamics in physics are 
based on assumption of a total energy conservation. Even more, the word "Hamiltonian" 
usually means "conservative" , while Lagrangian formalism is usually proved to be equivalent 
to Hamiltonian (therefore also conservative), as derived from a conservative Lagrangian 
energy function (for a comprehensive review see, e.g. PQ). The straightforward application 
of Hamiltonian/Lagrangian formalisms to human biomechanics would naturally inherit this 
conservative assumption [f| 

And this works fine for most individual movement simulations and predictions, in which 
the total human energy dissipations are insignificant. However, if we analyze a lOOm-dash 
sprinting motion, which is in case of top athletes finished under 10 s, we can recognize a 
significant slow-down after about 70 m in all athletes - despite of their strong intention to 



1 Ricci flow is a current hot topic in pure mathematics for which the latest Fields Medal was awarded to 
G. Perelman for the proof of 100 year-old Poincare Conjecture. 

2 An engineering- type approach to this problem would be simply adding dissipation and forcing into 
equations of motion, without deriving them form some "dissipative Lagrangian" and/or "dissipative Hamil- 
tonian"; this yields an extended Hamiltonian and/or Lagrangian formalism (see [2113114115116117118119]). 
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finish and win the race, which is an obvious sign of the total energy dissipation. This can 
be seen, for example, in a current record-braking speeddistance curve of Usain Bolt, the 
world-record holder with 9.69 s 10J, or in a former record-braking speeddistance curve of 
Carl Lewis, the former world-record holder (and 9 time Olympic gold medalist) with 9.86 
s (see Figure 3.7 in [9]). In other words, the total mechanical + biochemical energy of a 
sprinter cannot be conserved even for 10 s. So, if we want to develop a realistic model of 
intensive human motion that is longer than 7-8 s (not to speak for instance of a 4 hour 
tennis match) , we necessarily need to use the more advanced formalism of time-dependent 
mechanics. 

Similarly, if we analyze individual movements of gymnasts or pirouettes in ice skating, 
we can clearly see that the high speed of these movements is based on quickly-varying 
mass- inertia distribution of various body segments (mostly arms and legs). As the total 
mass-inertia matrix of a biomechanical system corresponds to the Riemannian metric tensor 
of its configuration manifold, we can formulate this problem in terms of time-dependent 
Riemannian geometry [H [I] . 

The purpose of this paper is to introduce a general framework for time-dependent biome- 
chanics, consisting of: 

1. human biomechanical configuration manifold and its (co)tangent bundles; 

2. biomechanical jet spaces and prolongation of locomotion vector-fields developed on the 
biomechanical configuration manifold; and 

3. time-dependent Lagrangian dynamics using biomechanical jet spaces. 

In addition, we will show that Riemannian geometrical basis of this framework is defined 
by the Ricci flow. In particular, we will show that the exponential-like decay of total 
biomechanical energy (due to exhaustion of biochemical resources [9]) is closely related to 
the Ricci flow on the configuration manifold of human motion. 

2 Human Biomechanical Manifold and its (Co)Tangent 
Bundles 

2.1 Humanoid Robot Dynamics 

Recall from [6] that representation of an ideal humanoid-robot motion is rigorously defined 
in terms of rotational constrained S'0(3)~groups in all main robot joints. Therefore, the 
configuration manifold Q ro b for humanoid dynamics is defined as a topological product of all 
included SO(3) groups, Q ro b = EL SO(3) 1 . Consequently, the natural stage for autonomous 
Lagrangian dynamics of robot motion is the tangent bundle TQ ro b, defined as follows. To 
each n— dimensional (nD) configuration manifold Q there is associated its 2nD velocity 
phase-space manifold, denoted by TQ and called the tangent bundle of Q. The original 
smooth manifold Q is called the base of TQ. There is an onto map ir : TQ — > Q, called 
the projection. Above each point x G Q there is a tangent space T X Q = 7r _1 (x) to Q at 
x, which is called a fibre. The fibre T. X Q C TQ is the subset of TQ, such that the total 
tangent bundle, TQ = [ | T X Q, is a disjoint union of tangent spaces T X Q to Q for all 

meQ 
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points x G Q. From dynamical perspective, the most important quantity in the tangent 
bundle concept is the smooth map v : Q — » TQ, which is an inverse to the projection n, i.e, 
ttov — Hq, tt(v(x)) = x. It is called the velocity vector-field. Its graph (x,v(x)) represents 
the cross-section of the tangent bundle TQ. This explains the dynamical term velocity 
phase-space, given to the tangent bundle TQ of the manifold Q. The tangent bundle is 
where tangent vectors live, and is itself a smooth manifold. Vector-fields are cross-sections 
of the tangent bundle. Robot's Lagrangian (energy function) is a natural energy function 
on the tangent bundle tqE 



2.2 Realistic Human Configuration Manifold 

On the other hand, human joints are more flexible than robot joints. Namely, every rota- 
tion in all synovial human joints is followed by the corresponding micro-translation, which 
occurs after the rotational amplitude is reached [6]. So, representation of human motion is 
rigorously defined in terms of Euclidean 5_E(3)-groups of full rigid-body motion [Til, 121 HI [T] 
in all main human joints (see Figure [1]). Therefore, the configuration manifold Q for hu- 
man dynamics is defined as a topological product of all included constrained SE(3) groups, 
Q = Yli SE(3y. Consequently, the natural stage for autonomous Lagrangian dynamics of 
human motion is the tangent bundle TQ (and for the corresponding autonomous Hamilto- 
nian dynamics is the cotangent bundle T*Q). 

Briefly, the Euclidean SE(3)-group is defined as a semidirect (noncommutative) product 
(denoted by >) of 3D rotations and 3D translations: SE(3) := 50(3) t> R 3 . Its most 
important subgroups are the following (for technical details see [H [T2l \T\): 



Subgroup 




Definition 




5*0(3), group of rotations 
in 3D (a spherical joint) 


Set of all proper orthogonal 
3 x 3 — rotational matrices 


SE{2) 1 special Euclidean group 
in 2D (all planar motions) 


Se 


t of all 3 x 3 — matric 
cos 9 sin 8 r x 
— sin 6 cos 9 r y 
1 


es: 


S0(2), group of rotations in 2D 
subgroup of 5£'(2)-group 
(a revolute joint) 


Set of all proper orthogonal 
2 x 2 — rotational matrices 
included in SE(2) — group 


R s , group of translations in 3D 
(all spatial displacements) 


Euclidean 3D vector space 



The configuration manifold Q = Y[i SE(3) 1 has the Riemannian geometry with the local 

3 The corresponding autonomous Hamiltonian robot dynamics takes place in the cotangent bundle 
T*Q rob , defined as follows. A dual notion to the tangent space T m Q to a smooth manifold Q at a point m 
is its cotangent space T^Q at the same point m. Similarly to the tangent bundle, for a smooth manifold Q 
of dimension n, its cotangent bundle T*Q is the disjoint union of all its cotangent spaces T^Q at all points 
m G Q, i.e., T*Q = | | T^Q. Therefore, the cotangent bundle of an n— manifold Q is the vector bundle 

T*Q = (TQ)* , the (real) dual of the tangent bundle TQ. The cotangent bundle is where 1-forms live, and 
is itself a smooth manifold. Covector— fields (1— forms) are cross-sections of the cotangent bundle. Robot's 
Hamiltonian is a natural energy function on the cotangent bundle. 
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Figure 1: The configuration manifold Q of the human musculo-skeletal dynamics is defined 
as a topological product of constrained SE(3) groups acting in all major (synovial) human 
joints, Q = H-SEfif. The manifold 

metric form: 

(g) = ds 2 = gij dx l dx-' , (Einstein's summation convention is in use) 

where g%j{x) is the material metric tensor defined by the biomechanical system's mass- 
inertia matrix and dx % are differentials of the local joint coordinates x l on Q. Besides giving 
the local distances between the points on the manifold Q, the Riemannian metric form (g) 
defines the system's kinetic energy: 

T = —giji x? , 

giving the Lagrangian equations of the conservative skeleton motion with kinctic-minus- 
potential energy Lagrangian L = T — V, with the corresponding geodesic form [5] 

j t L ± ,-L x ,=Q or F ■ l^J-'.r" 0. (1) 

where subscripts denote partial derivatives, while Tj k are the Christoffel symbols of the 
affine Levi-Civita connection of the biomechanical manifold Q. 

This is the basic geometrical structure for autonomous Lagrangian biomechanics. In the 
next section will extend this basic structure to embrace the time-dependent biomechanics. 

3 Biomechanical Jet Spaces 

In general, tangent and cotangent bundles, TM and T*M, of a smooth manifold M, are 
special cases of a more general geometrical object called fibre bundle, where the word fiber V 
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bundle Y*XxV 




Figure 2: A sketch of a fibre bundle Y rs X x V as a generalization of a product space 
X x 1/; left - main components; right - a few details (see text for explanation). 

of a map it : Y — > X denotes the preimage tt~ 1 (x) of an element x ^ X . It is a space which 
locally looks like a product of two spaces (similarly as a manifold locally looks like Euclidean 
space), but may possess a different global structure. To get a visual intuition behind this 
fundamental geometrical concept, we can say that a fibre bundle Y is a homeomorphic 
generalization of a product space X x V (see Figure [2| , where X and V are called the base 
and the fibre, respectively, it : Y — > X is called the projection, Y x = k (x) denotes a fibre 
over a point a; of the base X, while the map / = ir^ 1 : X — > V defines the cross-section, 
producing the graph (x,f(x)) in the bundle Y (e.g., in case of a tangent bundle, / = x 
represents a velocity vector-field). 

The main reason why we need to study fibre bundles is that all dynamical objects (in- 
cluding vectors, tensors, differential forms and gauge potentials) are their cross-sections, 
representing generalizations of graphs of continuous functions. 

By extending this line of formal bundle thinking, we come to the concept of a jet man- 
ifold, which is based on the idea of higher-order tangency, or higher-order contact, at 
some designated point on a smooth manifold. Namely, a pair of smooth manifold maps, 
fx, fi : M — > N (see Figure[3]), are said to be k— tangent (or tangent of order fc, or have a 
fcth order contact) at a point x on a domain manifold M, denoted by f\ ~ f^, iff 

fiix) called — tangent, 
d x h{x), called I — tangent, 
dxxfiix), called 2 — tangent, 
etc. to the order k 

In this way defined k— tangency is an equivalence relation, i.e., 

fi ~ h =>■ fi ~ /lj /l ~ fi ~ /3 =*> /l ~ /3j /i ~ /i. 

Now, a k—jet (or, a jet o/ order fc), denoted by j£f, of a smooth map / : M — > JV at a 
point a; G M (see Figure 13]), is defined as an equivalence class of fc— tangent maps at x, 

j k J : M -> JV = {/':/' is A; - tangent to / at x} 

In the same way as in the case of a map / : M — * N, x f(x), also in the case of the 
fc— jet j£f : M — > N, x i— > /(a;), the point x in the domain M is called the source o{ j%f 
and the point /(x) in the codomain N is the target oij^f. 



AO) = 
d x fi{x) = 
d X xfi{x) = 



5 




Figure 3: An intuitive geometrical picture behind the k— jet concept, based on the idea of a 
higher-order tangency (or, higher-order contact). 



We choose local coordinates on M and N in the neighborhood of the points x and f(x), 
respectively. Then the k— jet j£f of any map close to /, at any point close to x, can be 
given by its Taylor-series expansion at x, with coefficients up to degree k. Therefore, in a 
fixed coordinate chart, the A; —jet is given by: 

j£f : M — ► N = { collection of Taylor coefficients up to degree k }. 

The set of all A; —jets from M to N is called the k—jet space J k (M,N). It has a natural 
smooth-manifold structure. Also, a map from a k—jet manifold J k (M,N) to a smooth 
manifold M or N is called a k— jet bundle. 

For example, consider a simple function f : X — >Y, x t— > y = f(x), mapping the A— axis 
into the Y — axis in R 2 . In this case, A is a domain and Y is a codomain. At a chosen point 
x G A we have: 
a 0— jet is a graph: (x, f(x)); 
a 1— jet is a triple: (x, f(x), f'{x)); 
a 2— jet is a quadruple: (x, f(x), f'(x), f"(x)), 
and so on, up to the order k (where f'(x) = etc). 

The set of all k— jets from j£f : A — > Y is called the A:— jet manifold J fc (A, Y). 

We now turn back into the field of time-dependent human biomechanics, where the 
fundamental geometrical construct is the configuration fibre bundle n : Q — > R. Given a 
configuration fibre bundle Q — > M over the time axis R, we say that the 1—jet manifold 
J : (R, Q) is the set of equivalence classes j}s of sections s l : R — > Q of the bundle Q — > R, 
which arc identified by their values s'(t), as well as by the values of their partial derivatives 
d t s l = d t s l {t) at time points tel. The 1—jet manifold J X (R, Q) is coordinated by (t, x l ,x l ), 
so the 1-jets are local coordinate maps 

jls:R^Q, t^(t,x\x l ) 

Similarly, the 2— jet manifold J 2 (R, Q) is the set of equivalence classes j^s of sections s l : 
R — > Q of the configuration bundle 7r : Q — > R, which are identified by their values s l (t), 
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as well as the values of their first and second partial derivatives, d t s % = d t s l (t) and d u s l = 
d t tS l (t), respectively, at time points t € R. The 2-jet manifold J 2 (R, Q) is coordinated by 
(t, x l , ±\ x l ), so the 2-jets are local coordinate maps 



4. 



This global geometrical structure of time-dependent biomechanics is depicted in Figure 




time axis t e R 



Figure 4: Hierarchical geometrical structure of time-dependent biomechanics. Here, gener- 
alized coordinates q 1 , velocities q 1 and accelerations q 1 replace the corresponding rotational 
and translational joint coordinates a; 1 , velocities x 1 and accelerations x 1 . 



In the framework of biomechanics, we consider a pair of maps /i, /a ■ R — > Q from the real 
line R, representing the time t— axis, into a smooth n— dimensional Riemannian configuration 
manifold Q from Figure [1] We say that the two maps fi — fi (t) and = f% (t) have the 
same k— jet j t fe / at a specified time instant to € R, iff: 

1. /i(<) = /2(i) at t e R; and also 

2. the first fc terms of their Taylor-series expansions around to € R are equal. 

The set of all A: -jets j t fe / : R -> Q is the fc— jet manifold J k (R, Q). In particular, J X (R, Q) ^ 
R x TQ. 
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3.1 Jet Prolongation of Locomotion Vector-fields 

Consider an arbitrary locomotion vector-field u that is on the biomechanical configuration 
bundle tt : Q — > R defined by 

u = «*<9 t + u l (t, q j )d t . 

The so-called jet prolongation of the locomotion vector-field u onto the extended velocity 
phase-space, tht is the 1— jet manifold J 1 (M,(5), reads 

J 1 u = u t d t + u i d i + d t u i d t i . 

For example, consider a one-dimensional motion of a point particle subject to friction. 
It is described by the dynamic equation 

q = —kq, ( with k > 0). 

This is the Lagrangian equation for the Lagrangian function 

L = — exp[kt]q 2 dt. 

Let us consider the vector field 

v = d t - -qd q . 

Its jet prolongation J x v reads 

J l v = d t - ^qd q - ^q t d\. 

Generalization of all above jet structures to the k—jet manifold J fc (R, Q) is obvious, 

j^s :R->Q, t^(t, x\x\x\ x\ ...,x {k)l ). 

For more technical details on jet manifolds and bundles with their applications to mechanics 
and physics, see [H QH [H [HI [22] ) . 

4 Lagrangian Form of Time-Dependent Biomechanics 

The general form of time-dependent Lagrangian biomechanics with time- dependent La- 
grangian function L(t;x l ;x z ) defined on the 1— jet manifold X — J 1 (K, Q) = R x TQ, with 
local canonical coordinates: (t;x b ;x l ) — (time, coordinates and velocities) in active human 
joints, can be formulated as [H [1] 

^L ± i-L x i=J r i (t,x,x), (i = l,...,n), (2) 

where the coordinate and velocity partial derivatives of the Lagrangian are respectively 
denoted by L x i and L ± i. The right-hand side terms J-i(t,x,x) of © denote any type 
of external torques and forces, including [21 [3]: (i) excitation and contraction dynamics of 
muscular-actuators; (ii) autogenetic (spinal) reflex force controllers; (iii) rotational dynamics 
of hybrid robot actuators; (iv) (nonlinear) dissipative joint torques and forces; and (v) 
external stochastic perturbation torques and forces. 
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4.1 Time-Dependent Riemannian Geometry and Ricci Flow 



In the geodesic framework ([I}, the (in)stability of the biomechanical joint and center-of- 
mass trajectories is the (in)stability of the geodesies, and it is completely determined by the 
curvature properties of the underlying manifold according to the Jacobi equation of geodesic 
deviation [H [T] 

ds 2 3km ds ds 

whose solution J, usually called Jacobi variation field, locally measures the distance between 
nearby geodesies; D /ds stands for the covariant derivative along a geodesic and R? j km are 
the components of the Riemann curvature tensor. 

On the other hand, the mass-inertia matrix of human body segments, defining the Rie- 
mannian metric tensor gij — gij(x), need not be time-constant, that is, in general we have 
9ij — 9ij( x ) — 9ij{ti x )- Majority of fast movements in gymnastics are based on fast-changing 
mass distribution. This time-dependent Riemannian geometry can be formalized in terms 
of the Ricci flow equation (or, the parabolic Einstein equation), introduced by R. Hamilton 
in 1982 [21], that is the nonlinear heat-like evolution equatioro 

dtg%j = -2-Rjj, (3) 

for a time-dependent Riemannian metric g = gijif) on a smooth n— manifold Q with the 
Ricci curvature tensor This equation roughly says that we can deform any metric on 
a 2-surface or n— manifold by the negative of its curvature; after normalization, the final 
state of such deformation will be a metric with constant curvature. The factor of 2 in 
p| is more-or-less arbitrary, but the negative sign is essential to insure a kind of global 
volume exponential decay^ since the Ricci flow equation ^ is a kind of nonlinear geometric 
generalization of the standard linear heat equation 

d t u = Au. (4) 

Like the heat equation (j4|), the Ricci flow equation ([3|) is well behaved in forward time and 
acts as a kind of smoothing operator (but is usually impossible to solve in backward time). 

4 The current hot topic in geometric topology is the Ricci flow, a Riemannian evolution machinery that 
recently allowed G. Perelman to prove the celebrated Poincare Conjecture, a century-old mathematics 
problem and win him the 2006 Fields Medal (which he declined in a public controversy). The Poincare 
Conjecture can roughly be put as a question: Is a closed 3-manifold Q topologically a sphere if every closed 
curve in Q can be shrunk continuously to a point? In other words, Poincare conjectured: A simply-connected 
compact 3-manifold is diffeomorphic to the 3-sphere S 3 . 

5 This complex geometric process is globally similar to a generic exponential decay ODE: 

x = -Xf(x), 

for a positive function /(x). We can get some insight into its solution from the simple exponential decay 
ODE, 

x = —Ax with the solution x(t) = xoe~ At , 

(where x = x(t) is the observed quantity with its initial value xq and A is a positive decay constant), as well 
as the corresponding rath order rate equation (where ra > 1 is an integer), 

x = — Ax n with the solution = + (n — 1) Xt. 

~7i — I ™ 77 — 1 v ' 
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If some parts of a solid object are hot and others are cold, then, under the heat equation, 
heat will flow from hot to cold, so that the object gradually attains a uniform temperature. 
To some extent the Ricci flow behaves similarly, so that the Ricci curvature 'tries' to become 
more uniform [29j . thus resembling a monotonic entropy growth^ dtS > 0, which is due to 
the positive definiteness of the metric g^ > 0, and naturally implying the arrow of time [T]. 

In a suitable local coordinate system, the Ricci flow equation on a biomechanical 
configuration manifold Q has a nonlinear heat-type form, as follows. At any time t, we 
can choose local harmonic coordinates so that the coordinate functions are locally defined 
harmonic functions in the metric git). Then the Ricci flow takes the general form [3T] 

d t gij = A Q gij + Gij(g,dg), where (5) 

is the Laplace-Beltrami operator of the configuration manifold Q and Gijig, dg) is a lower- 
order term quadratic in g and its first order partial derivatives dg. From the analysis of 
nonlinear heat PDEs, one obtains existence and uniqueness of forward-time solutions to the 
Ricci flow on some time interval, starting at any smooth initial metric go on Q. 

The exponentially-decaying geometrical diffusion is a formal description for pirouettes 
in ice skating and fast rotational movements in gymnastics. 

4.2 Topological Lagrangian Dynamics 

The biomechanical jet space X — J 1 (R,X) = K x TQ gives rise to the fundamental 
n—groupoid, or n— category II n (X) (see [H Q]). In II n (X), 0-cells are points in X; 1— 
cells are paths in X (i.e., parameterized smooth maps / : [0,1] — > X); 2-cells are smooth 
homotopies (denoted by ~) of paths relative to endpoints (i.e., parameterized smooth maps 
h : [0,1] x [0,1] — ► X); 3-cells are smooth homotopies of homotopies of paths in X (i.e., 
parameterized smooth maps j : [0, 1] x [0, 1] x [0, 1] — > X). Categorical composition is defined 
by pasting paths and homotopies. In this way, the following recursive Lagrangian homotopy 

6 Note that two different kinds of entropy functional have been introduced into the theory of the Ricci 
flow, both motivated by concepts of entropy in thermodynamics, statistical mechanics and information 
theory. One is Hamilton's entropy, the other is Perelman's entropy. While in Hamilton's entropy, the scalar 
curvature R of the metric gij is viewed as the leading quantity of the system and plays the role of a probability 
density, in Perelman's entropy the leading quantity describing the system is the metric gij itself. Hamilton 
established the monotonicity of his entropy along the volume- normalized Ricci flow on the 2— sphere S 2 |26) . 
Perelman established the monotonicity of his entropy along the Ricci flow in all dimensions |30l . 
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dynamics emerges on the configuration manifold X: 

— cell : #o • in the higher cells below: t, s £ [0, 1]; 

L 



1 — cell : Xq • 



• X\ L : x — xi £ X, 



L : [0,1] -> X, L : x xi, xi = L(x ), L(0) = x , L(l) = xi; 
e.g., linear path: L(t) = (1 — t) xq + t Xi; or 
Lagrangian L — dynamics with endpoint conditions (xq,x\) : 
d 



dt 




(i = 1, ...,n); 



2 — cell : x i 



ft : [0, 1] x [0, 1] — ► ^T, /i : i — > L 2 , i 2 = ^(L^aro)), 
h(a;o,0) = L^aro), M^o, 1) = i 2 (^o), h(0,t) = x , h{l,i) = x x 
e.g., linear homotopy: h{xQ,t) = (1 — t) L(x ) + t L 2 (x )-. or 
homotopy between two Lagrangian (L , L 2 ) — dynamics 
with the same endpoint conditions (xq,x\) : 



— etc. 



and 



dt 



-Li 



with 



x(0) = x , 



= x x ; 



5 Conclusion 

In this paper we have presented time-dependent generalization of an 'ordinary' autonomous 
human musculo-skeletal biomechanics. Firstly, we have defined the basic configuration man- 
ifold Q of human musculo-skeletal biomechanics as an anthropomorphic chain of constrained 
Euclidean motion groups SE(3). Secondly, we have extended this base manifold by the real 
time axis K into the biomechanical configuration bundle Q — ► R. The time-dependent La- 
grangian biomechanics is defined on Q — > R using the formalism of first and second order 
jet manifolds as well as jet prolongations. Then we moved to time-dependent Riemannian 
geometry, governed by the Ricci-flow diffusion and showed that the exponential-like decay 
of total biomechanical energy (due to exhaustion of biochemical resources) is closely related 
to the Ricci-flow based geometrical diffusion. 
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